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>■ . The X-ray light curves of many gamma-ray bursts (GRBs) observed by the 

^ . Swift X-Ray Telescope (XRT) have a very steep-decay component (tail) following 

the prompt gamma-rays in the early phase and have some erratic flares occurring 
at a time from ~ 10^ up to ~ 10^ seconds. Based on the assumption that these 
prompt emission tails and flares are of "internal" origin and that their decline 
behaviors are dominated by the curvature effect of the fireball, we present a self- 
consistency test for this scenario with a sample of 36 prompt-emission-tails/fiares 
in 22 GRB XRT light curves. The curvature effect suggests that the temporal 
decay slope of the late steep-decay part of the light curves is a = 2 + /3, where 
13 is the X-ray spectral index. We derive the zero time (to) for each steep decay 
component by fitting the light curves with the constraint of a = 2+/5. Our results 
show that the to's of the prompt emission tails and the tails of well-separated 
flares are usually at the rising segment of the last pulse of the prompt emission or 
the corresponding X-ray flare, being self-consistent with the expectation of the 
internal dissipation models for the prompt emission and X-ray flares. Our results 
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indicate that each X-ray flare forms a distinct new episode of central engine 
activity and the GRB central engine remains active after the prompt emission 
is over, sometimes up to ~ 1 day after the GRB trigger (e.g. GRB 050502B & 
GRB 050724). This challenges the conventional central engine models and calls 
for new ideas to re-start the central engine. We further show that the on-set time 
of the late central engine activity does not depend on the GRB duration. We also 
identify a minority group of GRBs whose combined BAT-XRT light curves are 
smoothly connected without an abrupt transition between the prompt emission 
and the afterglow. These GRBs may have an external origin for both the prompt 
emission and the afterglow. 

Subject headings: gamma-rays: bursts — method: statistics 

1. Introduction 

The successful launch and operation of the Swift mission (Gehrels et al. 2004) have led 
to several important discoveries (e.g. Tagliaferri et al. 2005; Burrows et al. 2005; Gehrels 
et al. 2005; Fox ct al. 2005; Barthelmy et al. 2005; Roming et al. 2006; Cusumano et al. 
2006a). Combined analyses of the early data from the Burst Alert Telescope (BAT) and 
the X-Ray Telescope (XRT) for a large sample of bursts (Nousck et al. 2006; O'Brien ct al. 
2006) reveal a canonical X-ray afterglow lightcurve characterized by 5 components (Zhang 
et al. 2006a, see also Nousek et al. 2006): a steep-decay component associated with the 
GRB prompt emission tail, a shallow decay component likely due to refreshing of the forward 
shock, a normal decay component, a possible steep decay component following a jet break, 
as well as one or more X-ray flares. These new data provide unprecedented information to 
unveil the nature of these mysterious explosions. 

One of the outstanding problems in the pre- Swift era concerned the emission site of the 
prompt 7-ray emission. It is generally believed that GRB prompt emission originates at 
a distance internal to the fireball deceleration radius. The most widely discussed scenario 
is the internal shock model (Rees & Meszaros 1994; Kobayashi et al. 1997; Daigne & 
Mochkovitch 1998), but magnetic dissipation at an internal radius is also possible (e.g. 
Usov 1992; Thompson 1994; Drcnkhahn & Spruit 2002; Giannios & Spruit 2006). The 
broad-band afterglows, on the other hand, are produced by the external shocks when the 
flreball is decelerated by the ambient medium (Meszaros & Rees 1997; Sari et al. 1998). 
This "internal-|-external shock" model suggests that the prompt emission and the afterglow 
involve two distinct processes at two different emission sites. Alternatively, it has been 
argued that both the GRB prompt emission and their afterglows are produced in external 
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shocks, provided that the immediate medium near the burster is clumpy enough (Dermer 
& Mitman 1999; 2003). The evidence collected in the pre- Swift era cannot conclusively 
differentiate between the internal and the external models (see Zhang & Mcsaros 2004 for 
a critical review on the successes and limitations of both models). It is one of the scientific 
goals of the Swift to pin down the emission site of GRB prompt emission using early afterglow 
data. 

The steep-decay component commonly existing in early X-ray afterglows (Tagliaferri 
et al. 2005; Barthclmy et al. 2005b; Cusumano et al. 2006b; Vaughan et al. 2006) has 
been generally interpreted as the tail of the prompt emission (Zhang ct al. 2006a; Nousek 
et al. 2006; Panaitcscu ct al. 2006; Yamazaki et al. 2006; Lazzati & Begelman 2006). This 
component strongly suggests that the prompt emission and the afterglow are two distinct 
components, which supports the internal origin of the prompt emission (Zhang et al. 2006a; 
cf. C. D. Dermer, 2006, in preparation). The distinct X-ray flares typically show rapid rise 
and fall, with the ratio of the variability time scale and the epoch of the flare typically much 
less than unity, i.e. 6t/t <^ 1 (Burrows et al. 2005; Falcone et al. 2006; Romano et al. 2006). 
Burrows et al. (2005) proposed in their discovery paper that the flares are also produced in 
internal shocks at later times, which requires reactivation of the GRB central engine. Zhang 
et al. (2006a) performed more detailed analysis of various possible scenarios and concluded 
that the late internal dissipation model is the correct interpretation of X-ray flares. Similar 
conclusions have been also drawn by loka et al. (2005), Fan & Wei (2005), Falcone et al. 
(2006) and Romano et al. (2006) (cf. Piro et al. 2005; C. D. Dermer 2006, in preparation). 

An important clue to diagnose the internal origin of the prompt emission and the X-ray 
flares is the steep decay components following the prompt emission and the flares. These 
mark the sudden cessations of the emission, and the rapid decays are due to the observer 
receiving the progressively delayed emission from higher latitudes - the so-called "curvature 
effect" (Kumar & Panaitescu 2000; Dermer 2004; Zhang et al. 2006a; Fan & Wei 2005; 
Panaitescu et al. 2006; Dyks et al. 2006; Wu et al. 2006). The internal emission could 
be either from conventional internal shocks (Rees & Meszaros 1994; Kobayashi et al. 1997; 
Daigne & Mochkovitch 1998) or from magnetic dissipation at an internal radius (e.g. Usov 
1992; Thompson 1994; Drenkhahn & Spruit 2002; Giannios & Spruit 2006). Usually the de- 
tected emission contains the contributions from many emission episodes (e.g. emission from 
many internal shocks or magnetic dissipation regions). Each emission episode is expected to 
be followed by a curvature effect tail after the cessation of the emission (e.g. shock crosses 
the shell or magnetic reconnection flnishes). For highly overlapping emission episodes (e.g. 
the prompt emission), a curvature tail is usually buried beneath the rising lightcurve of an- 
other episode, so that it is difficult to observe a clear curvature effect signature. Although 
extensive studies have been made of the curvature effect in the prompt emission phase (Fen- 
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imore et al. 1996; Sari & Piran 1997; Ryde & Petrosian 2002; Norris 2002; Kocevski et al. 
2003; Qin et al. 2004; Shen et al. 2005; Qin & Lu 2005), no conclusive evidence supporting 
this effect has been presented (e.g. Kocevski et al. 2003). 

A clear, testable prediction of the curvature effect when the viewing angle is larger 
than the l/F cone (where F is the bulk Lorentz factor of the emission region) is that the 
temporal decay index a should be related to the spectral index (5 by the expression (Kumar 
& Panaitescu 2000) 

a = 2 + /3, (1) 

where the convention Fy oc t~"i'~^ is adopted. This relation is valid under the following 
assumptions, i.e. (1) there is a sharp drop off in the injection of accelerating electrons at 
a certain radius, after which the electrons cool adiabatically; (2) no spectral break energy 
crosses the observational band during the epoch of decay; and (3) the observational band is 
above the cooling frequency. When these conditions are satisfied, such a relation is rather 
robust regardless of the fireball history (e.g. Zhang et al. 2006a; Fan & Wei 2005; Dyks et 
al. 2006; Wu et al. 2006), and depends only on the jet structure when the hue of sight is 
outside the bright jet core (Dyks et al. 2006). 

Several effects could lead to deviations from eq.(l). The most important one is the 
zero time effect (io-effect). The io-effect is crucial, since when considering the multiple 
reactivation of the central engine, the GRB trigger time is no longer special (Zhang et al. 
2006a). Every time when the central engine is re-started, the new central engine time should 
be re-set as to- The conventional lightcurves are plotted as logF,^ - log(t — ^trigger), where 
^trigger IS the time when the GRB triggers BAT. When we consider the emission from a 
late central engine activity episode, the relevant decay slope (in order to be compared with 
the theoretical prediction in cq.[l]) is d\n Fy / dhi{t — tg) rather than din Ft, / d\n{t — ttrigger)- 
Properly shifting is therefore crucial to understand the real temporal decay index in 
the lightcurves^. For internal dissipation models (both internal shocks or internal magnetic 
dissipation) , the expected to of a certain flare should be at the beginning of the corresponding 
emission episode, i.e. near the starting point of the rising segment of the flare. Since both 
a and (3 could be directly derived from the Swift XRT observations, this provides a solid 
self-consistency test to the curvature-effect interpretation and the internal-origin hypothesis. 

Another effect that leads to deviations from eq.(l) is the overlapping effect (Zhang et 
al. 2006a), i.e. there is an underlying forward shock component beneath the steep decays. 
In order to perform a clean test to the curvature effect, such a component needs to be 
subtracted. Kumar & Panaitescu (2000) suggested that the contribution of the external 



The to issue is also relevant when discussing orphan afterglows from dirty fireballs (Huang et al. 2002). 
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shock component is significant if the ambient density is not lower than 10 ^ cm ^. 

In this paper we use the XRT data to test the curvature-effect interpretation and 
internal-origin hypothesis of the GRB prompt emission and X-ray flares. We take the GRB 
sample presented by O'Brien et al. (2006) and focus on the steep decay components fol- 
lowing the prompt emission and the X-ray flares. We assume that the curvature effect is 
the cause of the steep decays so that eq.(l) is valid ubiquitously. After subtracting the 
underlying contribution from the forward shock, we derive the values by fitting the light 
curves with our model. The object is to check whether the location of to is consistent with 
the expectation of the model, i.e. near the starting point of the relevant emission episode. 
The fitting model is described in §2. The data and the fitting results are presented in §3. 
Conclusions are drawn in §4 with some discussion. 



The rapid decay component is our primary interest. As discussed above we assume 
that this component is mainly contributed by the curvature effect. On the other hand, the 
overlapped contribution from an external shock component for n > 10^^ cm^"^ could result in 
a decay that would deviate significantly from eq.(l) (Kumar & Panaitescu 2000). Considering 
both the to-eSect and the overlapping effect, we model any steep decay component and the 
succeeding afterglow component with the function (Zhang et al. 2005) 



where f3 is the X-ray spectral index during the decay, to is the time zero point of the emission 
episode related to the decay (which in principle should be at the beginning of the rising 
segment of the last pulse of the prompt emission or the relevant X-ray flare if the curvature 
interpretation is correct), A and B are normalization parameters for both the rapid decay 
component and the underlying forward shock component, respectively, and C is the temporal 
index of the forward shock emission component. Please note that we take the zero point time 
of the external shock component (afterglow) as the GRB trigger time. This has been proved 
by detailed numerical simulations (Lazzati & Begelman 2006; S. Kobayashi et al. 2006, in 
preparation) . 



2. 



Model 
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3. Data and fitting results 

The Swift GRB sample we use is the same sample presented by O'Brien et al. (2006). 
These are the bursts detected by the Swift prior to October 1, 2005 for which prompt slews 
within 10 minutes were performed. The BAT data were processed using the standard BAT 
analysis software {Swift software v. 2.0). The XRT data were processed using xrtpipeline 
vO.8.8. The joined BAT-XRT hghtcurves were derived through extrapolating the BAT 
hghtcurves into the XRT band (0.3-10 keV) using the joint BAT-XRT spectral parameters. 
For the details of the data reduction procedure we refer to O'Brien et al. (2006). For our 
purpose, we include only those GRBs whose XRT hghtcurves have a steep decay component 
connecting to the prompt emission and those GRBs that harbor X-ray flares. In some bursts 
there are multiple flares, and hence, multiple steep decays. For these bursts we treat each 
steep decay independently. We consider only those prompt emission tails or flare tails whose 
decay slopes are steeper than -2 with the zero time set to the GRB trigger time. For the 
heavily-overlapped flares in the early XRT hght curves, we take only those that are well 
identified visually and without significant substructures. Our sample includes 36 prompt- 
emission/fiare tails from 22 GRBs, which are tabulated in Table 1. Their hghtcurves are 
collected in Fig.l. 

Technically one needs to identify the time interval for which our above test is performed. 
Since the steep-decay components (the tails of the prompt emission or the X-ray fiares) are 
required to satisfy the curvature effect, we need to search for a segment of lightcurve whose 
decay slope is steeper than -2. We choose a time window that contains 4 data points. We 
start the search from the beginning of the steep decay that follows the prompt emission or 
from the peak of a certain flare. We then move the time window to later times by shifting 
one data point in each step. We flt the decay slope for the four data points in the time 
window in each step until the slope becomes steeper than -2. By then the beginning of the 
time window is set to the starting point of the steep decay component. The end of the 
time interval for which our fit is performed is selected visually without a rigid criterion. For 
those prompt-emission/fiare tails with an underlying shallow decay (afterglow) component, 
we take the end of that component (before a further break if any) as the end of the time 
interval and flt the data using eq.(2). For some tails without superposition of an underlying 
afterglow component or otherwise highly overlapped with other flares, we simply choose the 
last data point of the steep decay component as the end of the time interval and flt the 
lightcurve using eq.(2) by setting B = 0. 

The GRBs 050406, 050502B, 050713B, and 050801 have a very flat segment following 
the tails of the prompt emission or flares. These flat segments have only a few data points, 
and we flx the C values in our fltting for these GRBs. Several bursts, e.g. GRBs 050712 (the 
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first and the second flares), 050713A, 050716 (the first fiare), 050730 (the 1st - 4th fiares), 
and 050822 (the first and second fiares) have heavily overlapped fiarcs before the typical long 
tail as observed in many other GRBs (e.g. 050126, 050219, 050315, etc.). We also fit these 
fiares. The XRT light curve of GRB 050714B has a rapid fiare in the tail segment. The level 
of this fiare is significantly higher than the tail. We subtract the fiare from the tail segment 
and fit the tail segment and the fiare independently. 

The spectral indices we use are from O'Brien et al. (2006). These spectral indices 
are derived from the spectral fitting to the overall X-ray data without considering spectral 
variability. In principle one should use the /? value during the steep decay to perform the 
test. However the photon counts during the steep decay only are usually too low to give a 
high significance fit. The X-ray spectra of some GRBs, such as GRB 050502B (Falcone et al. 
2006) and GRB 050607 (Pagani et al 2006), have detectable spectral variability. We check 
the difference made by this effect with GRB 050502B, and find that this effect does not 
significantly affect our fitting results. Throughout this analysis we use the spectral indices 
from O'Brien et al. (2006). 

Our fitting results are summarized in Table 1. In Figure 1, we mark the time interval 
of the fitting data and the fitting curve by sohd lines, and indicate the fitted to of each 
steep decay component by a vertical dash-dotted line. For those lightcurves with multiple 
tail/fiares, we identify the to of the i-th steep decay component as to,i- The reduced of 
the fits are also shown in Fig. 1. 

From Fig. 1 we find that, except those heavily-overlapped flares in the early XRT hght 
curves, the fltted to's of a good fraction of the well-identified tails are right at the beginning 
of the rising phase of the fiare or the last pulse of the prompt emission. These include GRBs 
050126, 050219A, 050319, 050406, 050422, 050502B (the late fiare), 050607, 050712 (the 
third fiare), 050716 (the second fiare), 050724 (the late fiare), 050803, and 050822 (the third 
fiare). This is well consistent with our starting hypothesis, i.e. the fiares are of internal origin 
and mark the re-activations of the central engine (Zhang et al. 2006a). We emphasize that 
the fitted ^o's are based on the hypothesis of the curvature effect as the origin of the steep 
decay, which is only relevant to internal models, not to external models. So even if some 
external models may also allow re-set of to before the fiares, the decay slope after the peak 
would follow some other predictions (typically fiatter than —(2 + (3), e.g. Wu et al. 2006). 
As a result, the fitted to's to satisfy those predictions would be significantly different from 
the ones we obtained, i.e. they should have a large off-set with respect to the rising part of 
the fiare. This means that one does not find a self-consistent solution for the external shock 
models. The impressive consistency displayed in the above bursts lends strong support to 
the internal origin of these X-ray flares. It is worth noticing that a very late flare (around 
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1 day) is evident in the XRT lightcurves of GRB 050502B (Falcone et al. 2006) and the 
short burst GRB 050724 (Barthelmy et al. 2005a). There have been questions as to whether 
the central engine can restart at such a late epoch and whether these features are due to a 
refreshed external shock origin (e.g. Panaitescu et al. 1998; Zhang & Meszaros 2002). Our 
results indicate remarkably that the to's for these two late flares are also right before the 
rising part of the flares. This gives strong support to the internal origin of these flares, and 
calls for central engine models that can operate as long as ~ 1 day. 

The fitted to's for some tails in GRBs 050315, 050712, 050713A, 050714B, 050716, 
050721, 050724, 050730, 050801, 050814, 050819 and 050822 arc all before the peak of the 
corresponding flare or the starting point of the steep decay. However, they arc not apparently 
located at the rising part of a flare or prompt emission pulse. These are all highly overlapping 
flares or continuously decreasing prompt emission, so that the rising segment of the flare or 
the prompt emission pulse is deeply buried beneath the continuous emission. These cases 
do not contradict the internal models and the curvature effect interpretations, although it 
does not directly support the scenario. Further complications may come from the possibility 
that the electron injection does not cease abruptly, that the observational band is below the 
cooling frequency, and that spectral breaks may cross the band during the decay, etc. In 
any case, it confirms that it is difficult to search for evidence of the curvature effect using 
the prompt GRB light curves. On the other hand, if one beheves the internal origin and 
curvature effect scenario (as is self-consistently tested in other bursts), then the fitted io's in 
these bursts give some indication of the central engine time when the flare or the last pulse 
of the prompt emission is powered. 

In Figure 2, we display the fltted to^s as a function of Tgo, the duration of the prompt 
emission. The cases of heavily overlapped flares are not included. It is found that the 
two quantities are not correlated, indicating that the on-set time of the late central engine 
activity does not depend on the duration of the prompt central engine activity. This suggests 
irregular, unpredictable behaviors of the GRB central engine. 

We have used the hghtcurve in the observer's frame to perform the above flts. Using the 
intrinsic lightcurves (i.e. systematically changing the time axis to the cosmic proper frame 
time t/{l + z)) shifts the searched time zero point to to/{l + z). So our conclusion does not 
depend on the time dilation effect since the lightcurves are simply compressed in the cosmic 
proper frame. 
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4. Conclusions and Discussion 

We have analyzed 36 prompt-emission/flare tails in the XRT light curves of 22 GRBs 
detected by the SwifthefoTe October 1, 2005 that show clear steep decay components. This 
is a sub-sample of that presented by O'Brien et al. (2006). Assuming that the tails are 
predominantly caused by the curvature effect, we derive the tg^s of these tails by fitting the 
XRT light curves with eq.(2). Our results (Fig.l and Table 1) suggest that usually the to^s 
are near the beginning of the rising segment of the last pulse of the prompt emission or the 
corresponding X-ray flare, which is consistent with the expectation of the internal dissipation 
models for the prompt emission and X-ray flares (e.g. Burrows et al. 2005; Zhang et al. 
2006a; Fan & Wei 2005; loka et al. 2005; Falcone et al. 2006; Romano et al. 2006). This 
suggests that the GRB central engine re-activates after the early prompt emission is over, 
sometimes up to days after the trigger (e.g. GRB 050502B & GRB 050724). This challenges 
the conventional central engine models and calls for new ideas to re-start the central engine 
(e.g. King et al. 2005; Perna et al. 2006; Fan et al. 2005; Proga & Zhang 2006; Dai et al. 
2006). We also show that the on-set times of the late central engine activity do not depend 
on Tgo, suggesting an erratic, unpredictable central engine. 

About one-third of the light curves in the O'Brien et al. (2006) sample do not show a 
distinct steep-decay component that connects the prompt emission and the afterglow. It is 
possible that in some GRBs (e.g. GRB 050525A, Blustin et al. 2006) the XRT observations 
started too late to catch the prompt emission tail. Nonetheless, a small fraction of XRT light 
curves show an apparently smooth transition from the prompt emission to the afterglow 
emission without a steep-decay bridge. The most prominent cases are GRB 050401 (De 
Pasquale et al. 2006), 050717 (Krimm et al. 2006), and 050826 (O'Brien et al. 2006). 
The joined BAT- XRT light curves of these GRBs (see O'Brien et al. 2006) have some 
low-amplitude flares overlapping on an otherwise smooth single power-law decay component 
extending to a late epoch. If these lightcurves are still interpreted within the frame work 
of internal-external model, as is favored by the majority of bursts in the sample (Fig.l), 
the prompt emission and the flares are of internal origin with a smaller emission efficiency 
than most other bursts (Zhang et al. 2006b). For example, the lightcurve of GRB 050717 
could be also well modelled by the superposition of a prompt emission tail emission and 
an underlying afterglow component (Krimm et al. 2006). Alternatively, it may be possible 
that the fireball in these cases is decelerated at a very early time so that both the prompt 
emission and the afterglow originate from the external shocks (e.g. Dermer & Mitman 1999; 
2003). This requires a high ambient density and a high initial Lorentz factor of the fireball. 
The low amplitude flares overlapping on the decaying lightcurves could be still from late 
internal dissipation of the central engine, or from some external shock related collisions (Wu 
et al. 2006), or else from possible density clumps in the medium (e.g. Dermer & Mitman 
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Table 1. The GRB sample and the fitting results 
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Table 1 — Continued 



GRB^ 


Interval(s)'' 
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*The number following the burst name denotes the sequence number of the steep decay components in the light curve according to the 
time sequence. 

''The time interval relative to the GRB trigger time for the data we use for the fitting, in units of seconds. 

"^Taken from O'Brien et al. 2006. These indices are derived from the spectral fitting to the overall X-ray data without considering the 
spectral variability. 

'^Thc decay part of the flare and the succeeding flat component have only four data points with large error bars. We fix the C value and 
fit the data. The error of to is thus small with a small reduced x^- 

''If the last data point at 704682 seconds is considered the reduced would be then ~ 1- This point drives most of the x^ values. 

^Some observational data points of this burst have an extremely small error, say, A log Fx < 0.1. The x^ of the fitting is mostly 
contributed by these data points, and the reduced x^ is unacceptable. To malje the fitting results more reasonable we talje A log Fx =0.1 
for those data points with A log Fx < 0.1. 

^These flares have few data points with large error bars. Their reduced x^ is too small. The fitting results are highly uncertain. 

''Excluding the superimposed X-ray fiare at 331-477 seconds. 

'Fit to the superimposed X-ray flare at 331-477 seconds only. 

jPit to the combined BAT-XRT data in 18-110 seconds. The flare at ~ 210 seconds does not satisfy our selection criteria (see in the 
text). 
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Fig. 1. — The combined BAT-XRT light curves with the best fitting results. The connected 
lines with error bars are the extrapolated BAT light curves. The solid circles with error bars 
represent the XRT observations. The thick solid curves are the best fits of our model, which 
cover the range of the fitted data. The vertical dotted lines mark the best fit to's- When 
more than one steep decay component is observed, the to of the i-th component is denoted 
as to^i. 




Fig 1. continued 
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Fig. 2. — The to as a function of Tgo for the tails/flares in our sample, except for those 
heavily overlapped flares in the early time XRT light curves. 



